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Abstract We study the friction coefficient of a macroscopic sphere in a viscous fluid at low Reynolds 
number. First, Kirkwood’s formula for the friction coefficient is reviewed on the basis of the Hamiltonian 
description of particle systems. According to this formula, the friction coefficient is expressed in terms 
of the stress correlation on the surface of the macroscopic sphere. Then, with the aid of large deviation 
theory, we relate the surface stress correlation to the stress correlation in the bulk of the fluid, where 
the latter is characterized by the viscosity in the Green-Kubo formula. By combining Kirkwood’s 
formula and the Green-Kubo formula in large deviation theory, we derive Stokes’ law without explicitly 
employing the hydrodynamic equations. 
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1 Introduction 

In equilibrium systems, macroscopic behavior is systematically described by a universal framework— 
thermodynamics [I]. On the basis of a microscopic description of equilibrium systems, the principle 
of equal a priori probabilities and Boltzmann’s principle reproduce the framework of thermodynam¬ 
ics, which is established as equilibrium statistical mechanics [2]. In contrast, there is still no theory 
describing the general behavior of non-equilibrium systems beyond the linear response regime. Thus, 
much effort has been devoted to the investigation of steady-state thermodynamics and non-equilibrium 
statistical mechanics [3J|31[51[n] • 

When liquids and gases are out of equilibrium but still remain in local equilibrium, their macro¬ 
scopic dynamical behavior is precisely described by the hydrodynamic equations [7]. A microscopic 
understanding of the hydrodynamic equations for the case of dilute gases was established through 
the Boltzmann equation [5], whereas it remains unclear for a general fluid. Non-trivial relations that 
are generally valid far from equilibrium, including the fluctuation theorems [iiiniiiiiiiiiiiiiiiiniiis] 
and the Jarzynski equality have been developed over the past two decades as a result of the 
time-reversal symmetry of microscopic mechanics. Thanks to such universal relations, we can easily re¬ 
derive certain well-known relations, such as the McLennan ensembles, the Green-Kubo formula, and 
the Kawasaki nonlinear response relation [mns]. Furthermore, by using a non-equilibrium identity 
similar to the fluctuation theorems and assuming a local Gibbs distribution at the initial time, the 
Navier-Stokes equation was derived for an isolated Hamiltonian system [^. Based on these achieve¬ 
ments, we believe this is an opportune moment to reconsider fluid dynamics from the viewpoint of 
statistical mechanics. 
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In this paper, we study the friction coefficient of a macroscopic sphere in a fluid. Starting from 
microscopic mechanics, Kirkwood first derived the linear response formula for the friction coefficient, 
7k, in the form |21j 

^ 3kBT Jo 

where fee is the Boltzmann constant, T is the temperature of the fluid, ( • )eq denotes a canonical 
ensemble average at temperature T, and Ft is the total force exerted on the sphere by the fluid at time 
t. The upper limit of the integral r should be much larger than the correlation time of the force Ft 
and much smaller than the relaxation time of the momentum of the sphere. This will be discussed in 
detail in Sect. |31 Note that r could go to infinity if we take the heavy mass limit for the sphere E1I131 
I24j . The theoretical results were confirmed by numerical experiments 

By focusing on a special class of fluids, we can obtain a more explicit expression for the friction 
coefficient. For instance, the friction coefficient for a dilute gas, 7 dg, has the form pmHn] 

7dg = ^p7^^\/27rmfcB'r, (2) 

where p is the number density of gas particles, TZ is the radius of the sphere, and m is the mass of 
a gas particle. As another typical example, the friction coefficient for a viscous fluid at low Reynolds 
number, 75 , is given by [7] 

7S = CpU, (3) 

where 77 is the viscosity of the fluid. The numerical coefficient C is equal to dvr for the slip boundary 
condition where the velocity of the fluid normal to the sphere at the boundary is equal to that of the 
sphere in this direction and the shear stresses on the sphere are equal to zero, and Gtt for the stick 
boundary condition where the velocity of the fluid at the boundary is equal to that of the sphere. The 
expression m is known as Stokes’ law. 

Kirkwood’s formula can be derived on the basis of a mechanical system, meaning that the validity 
of the formula is generally independent of the nature of the fluid within the linear response regime. 
Thus, Kirkwood’s formula ([T]) corresponds to (HI for a dilute gas, which was confirmed by Green [5T] . 
Furthermore, © corresponds to m for a viscous incompressible fluid at low Reynolds number, as 
confirmed by Zwanzig |32j on the basis of fluctuating hydrodynamics [71133] and Faxen’s law [33] in the 
case of the stick boundary condition. Note that the linear response formula is effective because the 
Reynolds number is low. 

We now consider the derivation of Stokes’ law m from Kirkwood’s formula o without explicitly 
employing the hydrodynamic equations. According to the Green-Kubo formula |35) . the viscosity of 
the fluid, 77 , can be expressed in terms of the stress correlation in the bulk of the fluid as 

j (4) 

where t) is the cc-component of the force on the unit area perpendicular to the y-axis at position 

r and time t (this is the xy component of the stress tensor and rg is an arbitrary position in the 
bulk of the fluid. Note that r is much larger than the correlation time of cr“^, and much smaller than 
the relaxation time of the momentum density field of the fluid. Precisely speaking, the condition of 
this r is different from that of r in Kirkwood’s formula, but we assume that there exists some r that 
satisfies both conditions. Because the right-hand side of ® can be expressed in terms of the stress 
correlation on the surface of the sphere, we can obtain Stokes’ law ([31 from Kirkwood’s formula (|T| and 
the Green-Kubo formula (|3| if the surface stress correlation is related to the bulk stress correlation. 

The main contribution of this paper is to establish the connection between the bulk and surface 
stress fluctuations. The basic concept is simple. The probability density of surface stress fluctuations 
is obtained from the probability density of bulk stress fluctuations by integrating out the other degrees 
of freedom. This procedure can be conducted in an elegant manner with the aid of large deviation 
theory. 

The remainder of this paper is organized as follows. In Sect. [3] we explain the setup of our model, 
and we review Kirkwood’s formula from the Hamiltonian description of particle systems in Sect. |3| 
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We derive the probability density of bulk stress fluctuations under two phenomenological assumptions 
in Sect. 01 We then express the probability density of surface stress fluctuations by integrating the 
probability density of bulk stress fluctuations. We apply a saddle-point method to this expression, 
and obtain the exact form of the probability density of surface stress fluctuations. By combining the 
obtained expression with Kirkwood’s formula, we derive Stokes’ law. These highlights are presented in 
Sect. [5l The final section is devoted to a brief summary and some concluding remarks. 

Throughout this paper, the superscripts a, b and a, /3 represent the indices in Cartesian coordinates 
(a;, y, z) and spherical coordinates (r, 0, ip), respectively, where {x, y, z) = (r sin 9 cos p, r sin 9 sin p, r cos 9) 
with r>0, O<0<7r, 0<<p< 27r. In addition, we employ Einstein’s summation convention for re¬ 
peated indices appearing in one term. 


2 Model 

We provide a three-dimensional mechanical description of our setup. The system consists of N bath 
particles of mass m and radius rbp, and one macroscopic sphere of mass M and radius 7^ in a cube of 
side length L. We assume that M and TZ are much larger than m and rbp, respectively. For simplicity, 
periodic boundary conditions are assumed. Let {ri,pi) {1 < i < N) be the position and momentum of 
the ith bath particle, and (i2, P) be those of the sphere. A collection of the positions and momenta 
of all particles is denoted by T = ... ,rN,pN,R,P), which represents the microscopic state of 

the system. 

The Hamiltonian of the system is given by 

fi |2 

= E +E - ^^-1)+- ^1) 

i—\ j>i 

where ^int is a short-range interaction potential between two bath particles, and ^sp is that between 
a bath particle and the sphere. We assume 


2M ’ 


(5) 


^sp(x)-too, asx-tT^ + rbp, ( 6 ) 

^sp(x)=0, for X > ^ + ?'bp+ Co, (7) 

where .^o — ''’bp- Pt denotes the solution of the Hamiltonian equations at time t for any state T at 
t = 0. In this setup, the energy is conserved, i.e., 


Hin) = H{r), 


and Liouville’s theorem 


m 

dr 


= 1 


holds. The total force acting on the sphere is given by 


F(T) 


dH{r) 

dR 

dPspilfi ~ -R|) 


( 8 ) 

(9) 


( 10 ) 


For convenience, we abbreviate J’(A) as Ft- The equation of motion for the sphere is written as 
dtPt=Ft. 

We assume that the system is initially in equilibrium at temperature T. Then, the initial probability 
density of F is given by 

Hjr) - Peg 
k^T 

where S'eq is the normalization constant. 


/eq(T) = exp 


( 11 ) 
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3 Kirkwood’s formula 

3.1 Derivation 

To consider the relaxation of the momentum of the sphere, we apply an impulsive force MVd{t) to the 
sphere when the system is in equilibrium at temperature T. Then, the probability density just after 
t = 0 is represented by 

H{r)- p -v -^{vy 


/o(T;y)=exp 


k^T 


( 12 ) 


where tf'(V') is the normalization constant and /o(T;0) = /eq(r). Using (|S]) and dH]), the probability 
density at time t is obtained as 


ft{r) = fyr.t-v) 

H{r) - p_,. y - ^{v) 


= exp 


kBT 


(13) 


In the following, the expectation values with respect to /o(T ; V) and ft{r) are denoted by ( • and 
( • )t, respectively. In particular, ( • )q corresponds to ( • )eq. 

We derive the exact formula for the friction coefficient for any V . Using (1131) and dtP-t = —F-t, 
we obtain 


Mn = foir-v)+ f dsdsfsin 
Jo 

= fyr-v)-j^dsfs{r) 


F-s V 


In addition, because F is independent of the momenta, we obtain 


Then, (fTTll and (fTKl) lead to 


= - 


kBT 


{F)o = iF)eo = 0 - 


J dr J^ds fsir)F<^F^_,V^ 


(14) 


(15) 


J dr_,/o(r_.;y)F“P^u' 


with 




itW) ^ 


1 


kBT . 


ds {f:f^)^ 


(16) 


(17) 


which is valid for any V. Note that is not equal to Ft /M in (I16L and thus the time dependence of 
(F“)t is described by that of yf^CU). Furthermore, in general, (P“)t is a nonlinear function of V. 

Next, we focus on the linear response regime. The dependence of 7 “^ on V is neglected. Within 
this regime, we can rewrite (ED as 

yf (0) = 7*5“' (18) 

with 


1 


It = 


kBT 7 
3fcBT, 


ds 


ds (F 7 F“), 


(19) 


where we have employed the isotropic property of the system in equilibrium. The expression m is 
the linear response formula for the friction coefficient, first derived by Kirkwood [21) . 
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and Tjr 


, respectively. 

V 


3.2 Time dependence of 

We denote the correlation time of Ft and the relaxation time of P by r„ 

We assume the separation of time scales represented by Tmicro ^ Tmacro- Here, (HSl) and {P)X = MV 
lead to 

= J dr fair-, V) (P“ - MV“) F’’ 

= 0 . 


( 20 ) 


Using dtPt = Ft and (EOl) , we can rewrite (inD as 

1 


iriv) = 


{PtFX - 


ksT 1 
ksT \ * ^0 


( 21 ) 


Thus, we find = 0 for any V when t :s> Tmacm, because the correlation between and F^ 

is considered to take about Tmacro- Furthermore, by considering (inD, we obtain 7 o^(V) = 0 for any 
V. Keeping these two limiting cases in mind, we conjecture the following time dependence of 7 “^(V). 
7 “^(V') increases when 0 < t ~ Tmicro- After that, 7 i“^(V) remains constant when Tmicro "C t Tmacro- 
Eventually, 7 i“^(V) gradually tends to zero when t ~ Tmacro- By recalling (ITTl) . we can express this 
behavior as 


71 a T7ib\ 


In the linear response regime, we denote the constant value of jt during the time interval Tmicro “C 
t -C Tinacro as 7 - Using this particular value of 7 , Tmacro can be expressed as M/ 7 . This 7 is the linear 
friction coefficient. 


>0, 

for 0 <C i — Tmicro; 


= 0, 

for Tmicro ^ Tmacro; 

(22) 

<0, 

for T ~ Tmacro- 



3.3 Other expressions 

We hereafter focus on the equilibrium case in which F is chosen according to the canonical ensemble 
/eq(T), and consider a finite time interval [ 0 ,t] that satisfies Tmicro t <C Tmacro- Because the motion 
of the sphere can be ignored up to r, we assume that the center of the sphere is fixed at the origin. 
Throughout this paper, for any physical quantity A{r), we define the time-averaged quantity by 


Air) = 


1 


Then, we can rewrite (HU) as 


dt A{rt). (23) 

In this derivation, it should be noted that (F/F/)eq is a function of |f — s|, and that (F’/F^)eq = 0 

when Tmicro ^ i Tmacro- 

We can also express the force Ft by the surface integration of a stress cr/p: 

F“(U) = J df2 nV/p^(7^, 12; U), (25) 

where 17 is a solid angle and n = (sin 0 cos v?, sin 0 sin (^, cos 0). By direct calculation, we find 

5<^>sp(|n-i^|)■ 


N 


yti.r-r) = -Y, 


dR^ 


[R^ -r\)D{r-n,R) 


(26) 
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with ^ 

D{r; n, R)= [ S {r - r, - {R - r,)0 ■ (27) 

Jo 

Indeed, by substituting (05)) into (1^51) . we obtain m, where we have used the divergence theorem and 

(i?^ — r^)d^D{r] ri, R) = S{r — Vi) — 5{r — R). (28) 

Furthermore, the stress tensor can be defined in spherical coordinates by 


r XX 
‘-^sp 

^sp 

‘-^sp 



^sp 

1 - 

b 

^sp 

(jyy 

‘-^sp 

ay^ 

^sp 

= Q 

9r 

^sp 

00 

^sp 

^sp 

zx 

_ sp 

a^y 

^sp 

sp. 


(j^r 

^sp 

^sp 

^sp 




with 


Q = 


sin 0 cos cos 0 cos — sin ip 
sin 0 sin ip cos 0 sin p cos p 

cos 0 — sin 0 0 


where denotes the transpose of Q. We then obtain 

F^[r) j dQ [cos6» <;(7^, Q- F) - sinfl c^f;(7^, Q- T)] . 


(29) 


(30) 


(31) 


Here, it is convenient to define the stress (T*(F) that represents the z-component of the force per 
unit area on the surface of the sphere. From (1311) and F^{r) = 47r77.^(T*(F), we have 


(^) = ^ / dn [cos0 <;(7^, 12; F) - sm0 af;(7^, 12; F)] . 


The formula 


) = 

is rewritten as 


7 = 


(47r7^2 




(32) 

(33) 


2kBT " ' 'sq 

The linear friction coefficient is expressed in terms of the correlation of the time-averaged stress at the 
surface, d*. 


4 Stress fluctuations in equilibrium viscous fluids 

Our goal is to derive Stokes’ law © from Kirkwood’s formula (13311 without explicitly employing the 
Stokes equations that are a linearized form of the Navier-Stokes equations in the low Reynolds number 
limit. Stokes’ law (|31) contains the viscosity of the fluid rj. Instead of introducing ry as a parameter of the 
Stokes equations, we define rj from the stress fluctuations in the bulk of the fluid by the Green-Kubo 
formula ©■ Thus, to derive Stokes’ law, we must relate the bulk stress correlation to the surface stress 
correlation. In this section, we derive the probability density of coarse-grained time-averaged stress 
fluctuations in equilibrium viscous fluids. 


4.1 The Green-Kubo formula 

We first consider the microscopic expression of the stress cr“*’(r; F) in the bulk of the equilibrium fluid. 
For the momentum density of the bath particles 


N 


n‘^{r;F) = J2Pidir-r^), 


^(r;T) is defined to satisfy the following continuity equation: 

dtR^r; Ft) = d" [a“'(r; Ft) + <(r; A)] ■ 


(34) 


( 35 ) 










Derivation of Stokes’ Law 


7 


Note that the stress a‘^^{r;r) corresponds to the force per unit area r)n^ exerted on a virtual 

surface described by the unit vector from the outer side, where the unit vector is perpendicular to 
the surface and directed toward the outer region. cr“p was given in (l26ll . Here, cr“^(r;H) is calculated 
as 

m 

i=l 

It can be directly confirmed that 

a^\r;r) = a^^{r-,r). (37) 

Note that, when there are interactions between bath particles across the periodic boundary, the right- 
hand side of (ESI) must be modified to take a non-zero value on the shortest line connecting two 
interacting bath particles. 

By employing the microscopic expression of the stress, the Green-Kubo formula (0]) is precisely 
expressed as 

= dt J d^r{a^y{ro-r)a^y{r;rt))^^, (38) 

where r is chosen such that r <C T^acro- T'm'j'cro is the correlation time of , and is 

the relaxation time of JT. In general, this r is different from that in Kirkwood’s formula (1331) . Here, 
we assume that the same value of r can be chosen in the two formulas, and use the same notation. 
One may notice the long-time tail of the stress correlation, which apparently breaks the separation of 
time scales. We assume that this long-time tail only appears in the regime t > because it originates 
from long wavelength fluctuations of locally conserved quantities. Therefore, we continue the argument 
without considering the long-time tail problem. 


d^int{\ri -rj\) 


d\ri 




D{r;r„rj). (36) 


4.2 Macroscopic fluctuation theory 


Let ^micro be the largest length scale appearing in the molecular description, and ^macro be the minimum 
length characterizing macroscopic behavior. In the case of viscous liquids, ^micro may be estimated as 
and ^macro 1® given by 77.. Because ^micro Cmacroi can choose H to satisfy t^micro ^ Cmacro- 
We then define the coarse-grained time-averaged stresses by 

Jv(r) 

with V(r) = [r“ — Al2,r°‘ + A/2]. In terms of the Green-Kubo formula (1551) is expressed 

as 

T 

We also note that 

(41) 

for \r — r'\ > A, because the correlation length of the stress field is on the order of ^micro- The relations 
dini) and m can be summarized as 


(d-^(r)d“^(r')),, = -^Si\r - r'|), (42) 

where 5((i) = 1 for d <g; H and S{d) ~ 0 for d> A. 

We now formulate a macroscopic fluctuation theory for the coarse-grained time-averaged stresses 
(t“^. Because d°'^(r) are obtained by integrating the microscopic stress over the region V(r), it is 
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reasonable, by considering the central limit theorem, to expect these to obey a Gaussian distri¬ 

bution. Here, we adapt the continuum description with the space mesh (ultraviolet cutoff) A. In this 
description, S{\r — r'\)/^ 6 {r — r'), and thus (1421) is rewritten as 

(d"^(r)d"^(r')) 3 q = - r'). (43) 

Next, we consider the statistical properties of the stresses. We start with the following decomposition 
of 0 -“^ into traceless and trace parts: 


,T“^(r) = -p{r)S^^ + s‘^\r) 


(44) 


with 


p(r) = 


a^^{r)+ayy{r)+a^^{r) 

3 


(45) 


where p and are called the mean pressure and the stress deviator tensor, respectively. Because 
is traceless, we have 

s^^r) + syy{r)+s^%r)=0. (46) 


Using the isotropic property and the assumption that the correlation length of the stress fluctuations 
is much less than H, we can express 


s^\r)r'^'{r')) = 

/ eq 


-b ^2(5“^'(5^“' -b 


ab ca'b' 


(5(i— r'), 


where the constants Bi, B 2 , and B 3 are determined below. By recalling (157)) . we obtain 
and then we find that 


which leads to Bi 


(r^(r)s^^(r')),q = (5^"(r)r^(r'))e,, 

= H 2 . Using (l46ll . we also have 

((s“(r) + syy{r) + s""(r)) (^"(r') + F^(r') + s""(r')))eq = 0, 


(47) 

(T^“(r), 

(48) 

(49) 


which gives B 3 


—2Bij3. Finally, the Green-Kubo formula (1551) leads to 

D 2fcB'r77 

i>i = -■ 

T 


(50) 


We summarize these relations as 


/ja&(^)^a'b'(^/)\ ^ 

\ / eq 


2 kBTr] 


A 


aba'b' 


S{r — r') 


(51) 


with 


^aba b' _ ^aa ^bb' _j_ ^ab' ^ba _ _^ab^ab' 


(52) 


Condition (1571) is necessary, but the statistical distribution of is not completely determined. 
The problem is to derive the statistical property of p. Because we are studying incompressible fluids, 
the pressure fluctuation arising from the density fluctuations need not be taken into account. Rather, 
we assume that p is determined from the balance of time-averaged forces in each region, which is 
expressed as 9^0-“*'(r) = 0. We express this in coordinate-free form as 


V • 


cr ir) 


= 0 . 


(53) 


It is obvious that the expectation values of the stresses satisfy the balance condition in the equilib¬ 
rium cases. However, assumption (1551) means that fluctuating stresses are already balanced in the 
macroscopic fluctuation theory. 
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4.3 Probability density of stresses 


For later convenience, we express ^ = Z\b with g{ab) = i and g{a'b') = j, where g{xx) = 1, 
giyy) = 2 , g{zz) = 3, g{xy) = g{yx) = 4, g{yz) = g[zy) = 5, and g{zx) = g{xz) = 6. That is, A is 
interpreted as a 6 x 6 matrix. (We use the same symbol A without confusion.) Explicitly, the matrix 
A is expressed as 


' 4/3 -2/3 -2/3 0 0 0 
-2/3 4/3 -2/3 0 0 0 
-2/3 -2/3 4/3 0 0 0 
0 0 0 1 0 0 
0 0 0 0 1 0 
0 0 0 0 0 1 


(54) 


As the probability density of is expressed using the inverse of A, we calculate the fcth eigenvalue 
Ck and corresponding eigenvector (j)]. defined by 




j=i 


(55) 


We find that 


{ 0 for fc = 1 
2 for fc = 2, 3 
1 for fc = 4, 5,6 


(56) 


The zero-eigenvector is calculated as (j){ = I/a/S for j = 1, 2,3 and (j){ = 0 for j = 4, 5,6. Although the 
matrix A is singular, we can define the pseudo-inverse of A as 


6 


^(A-i)*^</i=0, 

i=i 


(57) 


and 


6 


i=i 


for fc > 2. A straightforward calculation yields 


(58) 




' 1/3 -1/6 -1/6 0 0 0 
-1/6 1/3 -1/6 0 0 0 
-1/6 -1/6 1/3 0 0 0 
0 0 0 1 0 0 
0 0 0 0 1 0 
0 0 0 0 0 1 


(59) 


We identify (/\-i)“f'“ b with {A~^Y^ hy i = g{ab) and j = g{a'b')- 

Now, we write the probability density of {p(r), s“^(r)} such that the statistical properties of the 
macroscopic stress fields {p(t’), s“^(r)} are reproduced. Because the macroscopic stress fields obey the 
Gaussian distribution on the restricted configuration space given by (PI) and (ESI), the expression is 


V ({p(r),s“^(r)}) = Cexp 


j cPr s“''(r)(Z\-i)“'’“'^'s“''''(r) 


X n X n n ■ (^o) 


where C is the normalization constant. 
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Because the traceless condition (gSl) can be written as 


-ab. 


r) = 0, 


we find that b b jn dgO]) can be replaced by b ga b 

where ip'^ is an arbitrary vector. Explicitly, this is written as 


(61) 


(62) 


zi-i = Z\ 


-1 


form: 


2 qx 

9i + 92 

9i + 93 

94 

95 

96 

qi + <?2 

292 

92 + 93 

94 

95 

96 

<?1 + <?3 

92 + 93 

293 

94 

95 

96 

94 

94 

94 

0 

0 

0 

95 

95 

95 

0 

0 

0 

96 

96 

96 

0 

0 

0 

1/3, 

93 

= -1/6, 

and g 4 = 

95 

= 

76 = 


■ 1 

1/2 0 

0 

0 

0 





1/2 

1 0 

0 0 0 




1 _ 

0 

0 0 

0 0 0 





0 

0 0 

1 0 0 





0 

0 0 

0 1 0 





0 

0 0 

0 0 1 





(63) 


(64) 


Note that Zl ^ is the pseudo-inverse of Z\ given by 


Z\ = 


' 4/3 -2/3 0 0 0 0 

-2/3 4/3 0 0 0 0 

0 0 0 0 0 0 

0 0 0 1 0 0 

0 0 0 0 1 0 

0 0 0 0 0 1 


where A is obtained by setting the {zz) column and the {zz) row to zero in the matrix Z\ . 
The transformation of variables from to yields 


(65) 


V ({p(r'),(T“^(T')}) = C"exp [-Tl({p(r),(T“^(r)})] 

X ]^(5^i7'^^(r)-b (T^^(r)-b i7^^(r)-I-3p(r)^(5^V • d-(r)^ (66) 

r 

with 

I{{p{r),a‘^\r)}) = j d^r + p{r)5‘^^) (^- 1 )“''“'''' +p{r)5‘^’^') , ( 67 ) 

where C is the normalization constant. Because r >> Tmicro, 21({p(r), (7“^(r)}) corresponds to a large 
deviation function in probability theory [31137]. 

Finally, we consider the stress tensor in spherical coordinates. This is denoted by where a and 
/? represent r, d, or p. With a similar method, we obtain 

V ({p(r'), (T“'^(r)}) = C" exp [-tX ({p(r'), (T“'^(r)})] 

X d^d'’’’'(r) -b d^^{r) + + 3p{r)js(y ■ a (r)^ 


( 68 ) 
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with 


-I pOO pTT 

/ dr d0 dipr’^smO 

4/CBi rj jji 7q Jq 

X (,T“^(r) + p(r)5“^) (^-i)“/5“'/3' +p(r)5“'^') , (69) 

where C" is the normalization constant and we have taken the limit L —>■ oo. Z\ is obtained by setting 
the column and the row to zero in A . Note that the fourth-order isotropic tensor ^ 
takes the same form as in Cartesian coordinates. 

Because TZ ^ Cbp and the interaction potential between the sphere and a bath particle depends 
only on the distance between them, the momentum of the tangential direction at the surface of the 
sphere is conserved. This yields the boundary condition 

n) = n) = o. (70) 


5 Stress fluctuation at surface 


Let us return to Kirkwood’s formula (ESI). We want to calculate the friction constant 7 from the 
fluctuation intensity of d*. We first connect d* with 0 -“^. By recalling that TZ rbp and using dS]) and 
E3, we obtain 

= 0 (71) 

at the surface of the sphere. From © and EH), we obtain 

<p^(7^ + rbp+eo,f2;/^)=0. (72) 

The continuity of the total coarse-grained time-averaged stresses leads to 

d“'(7^, ^2) = 4 / df2' d“p'(7^, C'), (73) 

where S{TZ, f2) is the region on the surface of the sphere represented by 


27Z’^~^ 2TZ 


yl 


A 


tp- 


27^sin6l’‘^ ' 27^sin6» 


(74) 


Therefore, from (IS^ and (I7S1) . we find that cr* in (ISSl) is connected to by the relation 


(7 5 (t 


1 

dir 


dfi [cos 6 l d^’'(7^, n) - sin 6 » d•’'®(7^, C)] . 


(75) 


The probability density of cr* is determined from the statistical properties of in the bulk under 
conditions dZH and EH- This relation is formally written as 


P(d,) = [ VpVa'^f^ V ({p(r), d“'3(r)}) , 

J cr* :fix 


(76) 


where “cr*:fix” in the integral represents the condition given in (I75|) . 
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5.1 Saddle-point method 


We now evaluate the right-hand side of (1751) . In the asymptotic regime r ^ Tmicro, the functional inte¬ 
gral may be accurately evaluated by the saddle-point method. By introducing the Lagrange multiplier 
held A(r) = (A’'(r), A®(r), A‘^(r)) to take constraint (1551) into account, we obtain 


with 


and 


7((t*) = min 


P((T*) = C" exp [—r/(cr,)] 

1 


dr de 

^kB,Tr] Jq 

<7* :fix 


dip C ({d“'^(r)}) , 


C{{cr°‘^{r)}) =r2sin6» -(-A(r) • (v-V(r)) 


(77) 

(78) 

(79) 


where C" is the normalization constant, and s°'^{r) = a°‘^{r) + p{r)S°'^ with p{r) = + 

< 7 ^^{r) (T‘^‘^(r))/3. Note that the saddle-point method for the large deviation function has been 

rigorously verihed under certain conditions. This is called the contraction principle m in probability 
theory. In [Appendix A[ we present a derivation of the probability density of ct* within the framework 
of fluctuating hydrodynamics. 


^ _ y ^ / / / / 

As a reference for the argument below, we explicitly write V • <t and P gO‘ P ^3 


/ 

I dr 


r9 


V - (T = 


1 . da 
r de 


1 da'^^ 

r sin 0 d<^ 


-a 




d£^ 

dr 


1 da^^ _ 

r de rsin0 d^^ 


1 __|_ aO<^-ar<p 


r tan e 


V 


I 1 


1 darr , 2a^r 




+ ^ y 


r tan 6 

rO 


(80) 


dr 


r de 


r sin e dtp 


r tan e 


and 


+ {a^^f + {a^^f - - d^'^d’’’' 


+ (d’-'^)" + (d'^*^)" + (d‘^’-) , 

respectively, where we have used O, (|45|), and (El- 

Next, we calculate the right-hand side of (1751) . We start with the variation 


(81) 


poo pTT p2TT poo pTT p2'K p 

S dr d 6 dp C [{a°‘^{r)}) = dr dO d(^ < 9“ 
Jn Jo Jo Jn Jo Jo I 


dC 


d[d^'cF^P) 


5a°^^ 


dC 


da°‘P 


-d° 


dC 


a(i9“'d“^) 


f)r 

+ -^5X^ 
i9A“ 


The Hrst term of the right-hand side corresponds to the surface contribution. We find that 


/■2Tr 


dp d'^ 


because of the periodic boundary of p, and 


dC 


d(^Qva°‘P) 


,5d“^ 


[ de 


Jo 

s 

<b 


= 0 


= 0 , 


(82) 


(83) 


(84) 


dC 


d{d^a^P) 


= 0 


because 


(85) 
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at 0 = 0 and tt. Furthermore, as a result of cni), d’’® and are fixed to zero at the surface. Then, 

= 0 

at r = 7^. Thus, the first term of the right-hand side of (15^ becomes 

dL 


/ dr dO dip 

In Jo Jo 


d 


Sa^ 


Here, we impose 


dC 


d {d^d^^) 


= 0 


( 86 ) 

(87) 

( 88 ) 


r—lZ 


as the boundary condition of the variational problem. This is called a natural boundary condition, and 
explicitly gives 

X^(n,n)=0. (89) 

With this setup, the surface contribution of the variation vanishes, and we obtain the Euler-Lagrange 
equations 


dC 

dd°^P 


-d° 


dC 


a(9“'<T“/5) 


= 0 . 


(90) 


Using GSI), o, and dHU, we can write o as 

dX^ 

'19A® A' 


d’'’' = -p + 


d^^ = -p- 


d^v = _p. 


r do r 

1 dX'^ A’' 

r sin 6 dp r 


X^ 


rtan0 


^re 




d^r- ^ _ 


15A’ 

' dX^ 

X^\ 

r do 

dr 

r ) 

' 1 

dX^ 

laA^ 

r sin 

0 dp 

r do 

' 1 

dX^ 

dXv 

r sin 

0 dp 

dr 


(91) 




Note that (HSl) is equivalent to 

V • A = 0 (92) 

in the expression of (IM1) . By solving (IKHll . diTT]) . and (EH) with boundary conditions dZni), (ESI), and 
dMl), we obtain the probability density of cr*. The set of equations dSSl), Ell), and Ea coincide with 
the Stokes equations 


■qS/^u = Vp, 
V • ii = 0, 


(93) 

(94) 


when we set A(r) = 2qu{r). We may interpret u{r) as the macroscopic fluctuating velocity of the fluid 
generated by d*. 

By referring to the solution of the Stokes equations, we obtain 


A'-(r) 

A^(r) 


^1 — —^ 27?.ct* cosd, 

— 277.d* sind, 


A^(r) = 0, 

TZ^ _ 

Pir) = Poo - cosO, 


(95) 

(96) 

(97) 

(98) 
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where Poo is a constant. Then, (EH) leads to 


= -poo + 


cos 9 


a^\r) = a‘^‘^{r) = -p^, 
d’'«(r) = d®'^(r)=d^’'(r) = 0. 
Substituting these relations into (l78l) . we obtain 

(ct*)^ 27r7^^ 


2 A:bT?7 


This immediately yields 




k^Tp 


®i 2ttTZ^t 

Combining this with Kirkwood’s formula (IHHll . we arrive at Stokes’ law 


7 = 


(47r7^2)2 

2fcBr 




(99) 

( 100 ) 

( 101 ) 


( 102 ) 

(103) 

(104) 


5.2 Stokes’ law for a macroscopic sphere with a rough surface 


Because we have studied a macroscopic sphere with a smooth surface represented by a spherical 
symmetric interaction potential, the version of Stokes’ law we have obtained corresponds to that with 
the slip boundary condition in hydrodynamics. Indeed, we used the boundary condition (1701) . In reality, 
the surface of a sphere is assumed to be rough. Although the precise microscopic description of such a 
surface is not simple, it should be claimed that the boundary conditions dini) are not imposed in the 
macroscopic fluctuation theory. In this case, the natural boundary conditions 

A(7^,i7) = 0 (105) 


are imposed, so that the Euler-Lagrange equation can be obtained. 

By referring to the solution of the Stokes equations with these modified boundary conditions, we 
obtain 


3TZ TZ^\ 47?.cr*cos0 

A W=(l- 27+2-;)-3-. 

, / 3TZ 7?.^ \ 47?.ct* sin d 

AW = -(l-77-5;5j-3-' 

A*’(r)=0, 

_ 

Pir) =p'oo - cos6», 

where p'^ is a constant. Then, EH) leads to 



, /37^2 

27^4^ 

- I’oo + (^ ^2 

^4 J^+COS0, 

d^^(r) 

= ^^^{r) = -p'^ 

n* 

H-T-cr* cost*. 

d^^(r) 

= -t-ct* smt^, 

J.4 


d®^(r) 

II 

•i 

II 

O 



(106) 

(107) 

(108) 

(109) 

( 110 ) 
( 111 ) 

( 112 ) 

( 113 ) 
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Thus, we have that 




dfi 


Ak^Tr] 


7^2 7^4 

-V cr* COS U 


+ 


^ - • a 

—;-(T* siny 
r4 


= (o'*) 


27r 


4fcBT?7 Jt?, 
47r7^^ 

2 SfeeTr?’ 


dr r‘‘ 


27^4 47 ^® 107 ^®^ 

j,4 ^6 3 j,8 ) 


which leads to 




SkBTrj 


®'i AttTZ^t 

The substitution of (11151) into Kirkwood’s formula (1551) yields 

(47r7?.^)^r , 2\ ^ ^ 

^ = 2 tBT O®'* = 

This is Stokes’ law for cases with the stick boundary condition in hydrodynamics. 


(114) 


(115) 


(116) 


6 Concluding Remarks 

In this paper, we have derived Stokes’ law from Kirkwood’s formula and the Green-Kubo formula in 
the linear-response regime where the Stokes equations are valid. In this derivation, we did not assume 
the Stokes equations to describe Stokes flow, but rather formulated the relation between the stress 
fluctuations in the bulk and the stress fluctuation at the surface with the aid of large deviation theory. 

The heart of this derivation is the contraction principle, which is a standard technique in large 
deviation theory m- In the present case, we first used a phenomenological basis to apply the large 
deviation function to stress fluctuations in the bulk, and then applied the contraction principle. See [38] 
for another application of the contraction principle to the calculation of a large deviation function for 
time-averaged quantities. A similar argument has been employed to derive a generating function for the 
current fluctuation in stochastic non-equilibrium lattice gases [39) . Such a phenomenological argument 
is called the additivity principle, and the condition of its validity can be discussed within the framework 
of fluctuating hydrodynamics |40j . Similarly, our derivation can also be formulated within fluctuating 
hydrodynamics, as briefly explained in [Appendix A[ Although there is no one-to-one correspondence 
between our argument and that in the additivity principle, there may be a universal concept behind 
these two arguments. Determining a theoretical framework that provides a microscopic understanding 
of the phenomenological arguments would be an interesting problem. 

In this paper, we studied the force from a viscous fluid. More complex cases can be discussed 
using a similar formulation. For example, a cross-over from the dilute case to the viscous case 
which has been observed in numerical experiments [29] . may be one of our next targets. A simple 
but less trivial example may be a Brownian particle under a temperature gradient [51115511551155] . To 
elucidate the mechanism of the driving force and the friction force in this system, we must consider 
macroscopic fluctuation theory from microscopic mechanics, where the fluctuation of the energy flux 
should be taken into account [5^. Furthermore, the study of the force to a small system from an 
active environment is a hot topic in recent non-equilibrium statistical mechanics [46] . Examples of 
active environments include cytoskeleton networks (under chemical reaction) [57115511551150] . granular 
materials [5n[52l[53l[54] . and the assembly of small biological elements [551I561I571I58] . The law of the 
force from an active environment may be given by a phenomenological description in accordance with 
experimental observations. Flowever, the nature of the force is highly non-trivial. Thus, we believe that 
developing the theory by which the nature of the force from an active environment can be described 
would be a significant achievement. 

Finally, let us return to the formula described in (1751) . This claims that the large deviation function 
on the surface is determined from the variational principle associated with the large deviation function 
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in the bulk. Note that the large deviation function is called entropy, because the large deviation of 
the fluctuations of thermodynamic variables is equal to the thermodynamic entropy. Here, one may 
recall the so-called bulk-boundary correspondence in quantum many-body systems |591I601[611I621I63] . 
For the moment, we do not have any evidence for a direct connection with such theories, but it may 
be interesting to imagine such a possibility. 
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Appendix A: Fluctuating hydrodynamics 


We describe a viscous incompressible fluid at low Reynolds number by the fluctuating hydrodynamic 
equations in spherical coordinates. We denote by u{r,t), p{r,t), and the velocity, 

pressure, stress tensor, and random stress tensor at position r and time t, respectively. In this case, 
each component of the stress tensor is written as 


rr 

(7 




a 


dip 


a 


ipr 


-P' 

-P' 


du^ 


2^7 + ZT 


-p 2r] 


1 du^ . W 
r do r 
1 du'^ 
r sin 6 dp 
1 du"" du^ 

r do dr r 

1 du^ 1 du^f 
r do 
du^ 
dr 


„S9 


U 


6 


U 

-h 

r r tan 0 


„re 


r sin 0 dp 
1 du"" 
r sin 0 dp 


r tan0 


-I- s 


dip 


+ s 


ipT 


(A.l) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 


We assume that the temperature T and density p of the fluid are constant and homogeneous. Because 
we are focusing on incompressible fluid, we assume that the bulk viscosity C is equal to zero and that 
V ■ u = 0. This is written as 

du"" 2u’" Idu^ u® 1 du‘9 ^ 

-1-1-^ - 1 - n 4-(^-7) 

dr r r dO r tan 0 r sm 0 dp 

Furthermore, the time evolution equation of u{r,t) is assumed to be given by 


'’a "■ 


(A.8) 


because the Reynolds number is sufficiently low. Finally, the random stresses are assumed to be zero- 
mean Gaussian white noises with covariance 


where we have used ^ 


= 2fcBr?7A“^“'^' ~ ~ - t'), 


^apa'P' _ gaa gaP' ^pa _ _gapga'P' 

3 

0. Note that (IA.9|1 leads to 

s’'’' -k s^® -k s'^'^ = 0. 


(A.9) 

(A.IO) 


(All) 
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By recalling (TOI) . (TOj) . (TOl) . and (E21), we also have 

+ 3p = 0. (A.12) 

For any time-dependent quantity A(t), we denote its path during the time interval [0,t] by [A]. 
Then, using (IA.9I) and the Gaussian property of s“^, we obtain the probability density of {[p], 
in the form 


V {{[p], K]}) = Co exp [-rX ({[p], [u“]})] 

t r 

X n n ^ (a’'’' + -k 3p) (A.13) 

t r 

with 

l{{[p],K%[u‘^]}) = drr^J (A.14) 

where Cq is the normalization constant and in the right-hand side is related to (p, a°‘^, u°‘) through 
(lA.ll) (IA.6I) . By applying the contraction principle to (IA.13I1 . we obtain the probability density of the 
surface stress fluctuations as 


P(d,) = f VpVa‘^PVu‘^ V ({[p], K]}) 

J (7* :fix 

= Cpexp [-t/(o-*)] 


with 


-I pT poo pTT plTT 

I{^*)= min / dt dr d0 dp £ [m“]}) , 

{[ct“'’].[u“]} 4 A:bJ rjT Jq Jq Jq 


cr* :fix 


and 


C ({K^], K]}) = sine L«/3(^-i)“/3“'/3'go'/3' _ y • ^) -p As (V • It) 


(A.15) 

(A.16) 


(A.17) 


where Cq is the normalization constant, Ai(r,t) and X 2 (r,t) are the Lagrange multiplier fields, p = 
— (cr’'’' -I- cr®® -I- cr'^‘^)/3, and “CT*:fix” represents the condition given in 


1 

4ttt 


'0 


dt J df2 [ 


If a time-independent configuration with u 
as in ((77l) -d7^. 


cosd cr’'’'(7?., C, t) — sind cr’'®(77., 17, t)] . (A.18) 

= 0 is the minimizer of (IA.16|1 . we obtain the same result 
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